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We investigate thermoelectric efficiency of systems with broken time reversal symmetry under a 
three-terminal transport. Using a model of Aharonov-Bohm interferometer formed with three non- 
interacting quantum dots, we show that Carnot efficiency r\c can be achieved when the thermopower 
is a symmetric function of the applied magnetic field. On the other hand, the maximal value of 
the efficiency at maximum power is obtained for asymmetric thermopower. Indeed, we show that 
Curzon-Ahlborn limit is exceeded within the linear response regime in our model. Moreover, we 
investigate thermoelectric efficiency for random Hamiltonians drawn from the Gaussian Unitary En- 
semble and for a more abstract transmission model. In this latter model we find that the efficiency 
is improved using sharp energy-dependent transmission functions. 

PACS numbers: 72.20.Pa, 05.70.Ln 
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I. INTRODUCTION 

Thcrmoelectrics convert temperature gradients into 
electric voltages and vice versa. Strong demand for cost- 
effective pollution free form of energy conversion has re- 
sulted in a plethora of thermoelectric based applications. 
Unfortunately, the practical usage is limited by extremely 
low performance of thermoelectric materials. Increasing 
the efficiency of thermoelectric materials is one of the 
main lines of current thermoelectric research 1-7 . 

In the linear response regime, the performance of ther- 
moelectric materials is characterized by a single dimen- 
sionless parameter called figure of merit ZT, which is a 
combination of the main transport properties of a ma- 
terial, i.e. the electric conductivity a, the thermal con- 
ductivity k and the thermopower S, as well as of the ab- 
solute temperature T: ZT = (aS 2 /k)T. The maximum 
efficiency is given by 



Vc 



sJZT + 1-1 



VZT+T+1 



(1) 



Carnot efficiency r\c is reached in the limit of ZT — > oo. 
Although thermodynamics does not impose any upper 
bound on ZT , the strong interdependence between the 
electric and thermal transport properties makes it ex- 
tremely hard to increase the value of ZT above 1. To 
compete with the existing mechanical based energy con- 
version applications, development of thermoelectric ma- 
terials with ZT at least above 3 is required. An attractive 
alternative to increase the performance is to use broken 
time reversal symmetry systems for which the efficiency is 
determined by two parameters, an asymmetry parameter 
x and a "figure of merit" y generalizing ZT-. In these sys- 
tems, maximum efficiency r) max can be as high as Carnot 
efficiency even with very low values of figure of merit 



y, provided the asymmetry parameter x is very large. 
Note that in time reversal symmetric systems x = 1 and 
V ZT. 

In a thermodynamic system, Carnot efficiency is ob- 
tained for fully reversible transformations. This requires 
quasistatic transformations and consequently the derived 
output power is zero. Hence, the notion of efficiency at 
maximum power was introduced and in the linear re- 
sponse regime is given by^ 



?7(w max ) 



ZT 



2 ZT + 2 



(2) 



Note that in the limit of ZT —> oo, 77(w max ) takes the 
maximum value of r\c /2. This upper bound is commonly 
referred to as the Curzon-Ahlborn limit^r— . In princi- 
ple, for broken time reversal symmetry systems this limit 
can be exceeded with values of the asymmetry parame- 
ter such that |x| > 1—. Moreover, r/(u> max ) — ► rjc when 
|ar| — > oo (always within the linear response regime). 
Hence, it is potentially of practical relevance to study the 
thermoelectric transport in broken time reversal symme- 
try systems. 

Thermopower for a broken time reversal symmetry sys- 
tem is in general asymmetric with respect to the time re- 
versibility breaking parameter B i.e., 5(B) ^ S(— B). 
The asymmetry parameter x — S(B)/S(— B). In a 
non-interacting system, inelastic scattering can result in 
asymmetric thermopowe r 20 ' 21 . Conveniently, this can be 
achieved with the introduction of noise by means of a 
third terminal (probe) with its temperature and chemi- 
cal potential adjusted such that there is no net average 
flux of particles and heat between the terminal and the 
system. Various aspects of three-terminal thermoelec- 
tric transport were investigated in22r— Large asymme- 
try in the thermopower were obtained in Refj^ -, using 
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an Aharonov-Bohm interferometer model, first treated in 
the context of thermoelectric transport in Ref^i. How- 
ever, the obtained efficiency was very low^S. Also, follow- 
ing the above lines, a classical deterministic three termi- 
nal transport model was studied, showing large asymme- 
try with very low efficiency^. Hence, it remains to be 
seen whether and under what conditions efficiency close 
to Carnot can be obtained using broken time reversal 
symmetry system under a three terminal transport. 

In this paper, we investigate the thermoelectric effi- 
ciency of non-interacting systems with broken time re- 
versal symmetry under three terminal transport in the 
linear response regime. First we consider a model of 
Aharonov-Bohm interferometer formed with three dots 
as analyzed in Ref^S. In contrast to the earlier work, 
we focus on the efficiency obtained over the global opti- 
mization of all system and reservoir parameters following 
simulated annealing method^. Our results show that 
Carnot efficiency r\c can be obtained for maximum ef- 
ficiency r\raax when the thermopower is symmetric i.e., 
x = 1. By adding asymmetry, r\ max decreases. However, 
the efficiency at maximum power rj(uj max ) is maximum 
when there is asymmetry in thermopower (x ^ 1). To 
the best of our knowledge for the first time we show that 
the Curzon-Ahlborn limit r?(w max ) = ?7c/2, which is a 
rigorous upper bound for systems with time-reversal sym- 
metry, is exceeded in the linear response regime. How- 
ever, both the efficiency r\ max and ??(u; ma x) decrease with 
increase in asymmetry of thermopower at large values 
of |x|. Therefore, we study broader classes of models 
in an attempt to improve efficiency at large asymmetry. 
We consider random Hamiltonians drawn from Gaussian 
Unitary Ensemble (GUE) and finally an abstract model 
of transmission probabilities. While in this latter model 
we show that ry(w max ) can take values as high as Q.hlric, 
still we could not find, after global optimization of all 
parameters of the model, large efficiency at large asym- 
metry. Our results obtained for very broad classes of 
non-interacting models implies that it is practically very 
hard, if not impossible, to achieve in such models large 
efficiency with large asymmetry in thermopower. The 
reason of this still remains to be understood. 

The paper is structured as follows: In Sec. II, we re- 
view the calculations of efficiency of broken time reversal 
symmetry systems under three terminal transport. De- 
pendence of efficiency on asymmetry of thermopower for 
various models with broken time reversal symmetry is 
analyzed in Sec. III. Finally, Sec. IV summarizes our 
results. 



II. MODEL AND METHOD 

A. General setup 

The general set up consists of a system in contact 
with two reservoirs left (L) and right (R) at tempera- 
tures Tl = T + AT, Tr = T and chemical potentials 



fJ-L = /i + A/i, /ip = fi. Inelastic scattering effects are 
simulated by means of a third (probe) reservoir (P) at 
temperature Tp = T + ATp and chemical potential 
Hp = M + A/ip. Let J p k and J Eh denote the electric and 
energy currents from the fcth reservoir (k=L,R,P) into 
the system, with the steady-state constraints of charge 
and energy conservation J2k Jpk = Oj^fc ^Ek = 0. The 
sum of the entropy production rates at the reservoirs 
reads S = J2k(^ Ek ~~ HkJpk)/Tk- Within linear response, 
S = J • X = 5Zi=i Ji^-ii where J and X are four dimen- 
sional vectors defined as 



J — (eJpL, J q L,eJ p p, J q p), 
,A^ AT A^p ATp 
eT ' T 2 ' eT ' T 2 '' 



(3) 
(4) 



Here e is the electron charge and J q k = JEk — f^Jpk is 
the heat current. The relation between the fluxes Ji and 
the thermodynamic forces Xi within linear irreversible 
thermodynamics is 



LX. 



(5) 



where L is a 4 x 4 Onsager matrix, J and X are written 
as column vectors. Eq. ([5]) can be written in the block 
matrix form as 



J/3 



L/3ct 



X a 

X/3 



(G) 



where a stands for (1,2) and (i for (3,4). 

The probe reservoir is adjusted in such a way that 3p — 
J3 = J4 = 0, that is, the net electric and heat flow from 
the probe into the system vanishes. This implies that 
X/3 = — L^L/JqXq. and 



L'X n 



L' 



(7) 

Thus, the problem has been reduced to two coupled fluxes 
as 
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J 1 J 1 



X-2 



(8) 



where the reduced 2x2 Onsager matrix L' satisfies the 
Onsager-Casimir relations 



LUB)=L'(-B) (i,i = 1,2). 



(9) 



Here, B is a magnetic field breaking the time reversal 
symmetry. Seebeck and Peltier coefficients are given by 
S = L' u /(eTL' u ) and n = L' 21 /(eL' n ). Thermopower is 
asymmetric when L' 12 ^ L' 2 i, i.e., n ^ ST. 
Maximum efficiency is given by 



Vex 



(10) 



and depends on two parameters: the asymmetry param- 
eter x and the figure of merit y, where 



12 



T' ' 
^21 

^12-^21 
dctL' ' 



(11) 
(12) 
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Efficiency at maximum power is 

s Vc xy 



2 2 



y 



(13) 



B. Noninteracting systems 

Thermoelectric efficiency can be calculated exactly for 
noninteracting models by means of Landauer-Biittiker 
formalism^. Consider then a noninteracting system with 
Hamiltonian Hs- We model the reservoirs as ideal Fermi 
gases with distribution function fk(E) = {exp[(E — 
Hk)/kBTk\ + l} -1 - The electric and heat currents from 
the left reservoir are given by 

Jl = T [ dEj2[r k L(E)f L (E) -T Lk (E)fk(E)}, 

^ k 

J2 = \l dEj2(E-»L)[T kL (E)f L (E) 



-TLk(E)fk(E)], 



(14) 



where tm(E) is the transmission probability from reser- 
voir I to reservoir k at energy E. Analogous expressions 
can be written for J 3 and J4, provided the terminal L is 
substituted by P. 

The Onsager coefficients L{j are obtained from the lin- 
ear response expansion of the currents Ji as 



L n = — dE }2 T Lk(E) 

rp p + OO 

Li2 = — / dE(E- ^)J2r L k(E) 

n k+L 

rp p+OO 

L22 = T dE(E- ^) 2 J2r Lk (E) 

" J -00 ,,,, 



dfjE) 
dE 



k^L 



dfjE) 
dE 

df(E) 



dE 



L21 — L12 



(15) 



Analogous formulas are obtained for L33, L34 = L43, and 
£44, with the P terminal used instead of L. Similarly, the 
off diagonal block elements are obtained as 



L13 

L2i 

L23 
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/+00 
dEr LP (E) 
-OO 

/ + OO 
dE {E-^)t lp {E) 
-OO 

/ + OO 
dE {E - ii) 2 t lp {E) 
-OO 



df(E) 



df(E) 



dE 
df(E) 



dE 



(16) 



Using t pl (E) instead of t lp (E) in Eq. dTBJ), L 31 , L 32 = 
L41, and L42 are obtained. 

The transmission probabilities r pq (E) are calculated as 

Tpq (E) = Tr[T p (E)G(E)T q (E)G\E)l (17) 



where the broadening matrices Tk are defined in terms 
of the self-energies as Tk = i[£fc(-E) — S^E 1 )] and 
the (retarded) system Green function G(E) = [E — Hs — 
Ek^kiET 1 . ' 

Explicit expression for the matrix elements of the re- 
duced 2x2 

Onsager matrix L' can be derived from Eq. ([7]). We 
obtain 





= L n - 
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dI(E 3 
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dI{E 1 
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E X ){E 2 


-E 3 ). 



(18) 



Here dl and D arc 



dl = 



>. 2 T 3 
h 3 



U dE * 



.1=1 



df_ 

dEt 



x t lp {Ei ) [t pl (E 2 ) + t pr (E 2 )]t pl (E 3 ) , 
[r PL (Ei) + T PR (Ei)}} {Ex - ji){Ei - E 2 ). (19) 



III. RESULTS 

A. Aharonov-Bohm interferometer 

Here, we discuss the thermoelectric efficiency of a three 
quantum dot ring structure pierced by an Aharonov- 
Bohm flux and connected to three reservoirs, with each 
dot connected independently to one reservoir. The sys- 
tem is described by the Hamiltonian 



i0/3 



H S = ^ 6 k c l c k + (tLRC R C L e 



+ t RP c P c R e^ 3 + t LP c[c P e^ 3 + H.c), (20) 

where e k are the on-site energies, tij are the hopping 
strengths («,j = L, i?,P), <p is the flux. We assume the 
wide-band limit and hence the broadening matrices are 
given by Tk — ^k^Ck- Here, 7^ measures the tunneling 
rate of electrons between the reservoir k and the system. 

When there is anisotropy in the system (e/c ^ ey), and 
the Aharonov-Bohm flux <j) is non-zero, the off diagonal 
elements of the reduced Onsager matrix L' are asymmet- 
ric functions of the flux, i.e., L' 12 {4>) 7^ L' 21 {4>) = L' 12 {— 4>). 
The asymmetry parameter x defined in Eq. (Ill[) is the 
ratio of the off diagonal Onsager matrix elements and is 
in general different from the time reversal symmetric case 
where x = 1. 

For this model, the efficiency is a function of 12 in- 
dependent parameters: ei, e 2 , e 3l t LR , t RP , t LP , 4>, j L , 
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FIG. 2. Optimized values of efficiency at maximum power 
rj(u)max) versus asymmetry parameter x for the Aharonov- 
Bohm interferometer. The inset shows the variation of 
Ti(uJmax) near the symmetric value x = 1. Note that the 
Curzon-Ahlborn limit r]{ui m ax) = f?c/2 is exceeded in this 
regime. rj{ui max ) is scaled in units of r\c in this figure and 
thereafter. 



FIG. 1. Optimized values of maximum efficiency rjmax as 
a function of the asymmetry parameter x for the Aharonov- 
Bohm interferometer discussed in the text, with a linear [panel 
(a)] and a logarithmic [panel (b)] scale for J7 max . Note that 
rjmax is scaled in units of rjc in this figure and thereafter. 



7.R, IP, T, fi. We maximize the efficiency over all these 
parameters using the simulated annealing method. For 
numerical convenience, the parameters are restricted to 
the intervals e t € [-5,5], Uj € [-5,5], T G [1CT 3 ,100], 
H G [-5,5], <p G [—1,1], 7i G [10~ 3 ,1] (in units such 
that e = /i = fcs = l)- Note that the results are un- 
changed by varying the range of parameter values by a 
few times. Fig. []Ta) shows the dependence of the max- 
imum efficiency r\ max on the asymmetry parameter 
In the absence of asymmetry in the thermopower i.e., 
when x = 1, Carnot efficiency is reached. As the asym- 
metry is introduced, the optimized maximum efficiency 
is always less than Carnot efficiency. 

We also plot the dependence of the optimized maxi- 
mum efficiency r) max on a; in a lin-log scale in Fig. [IJb) . 
It is clear from the figure that for \x\ < 1, r\ max increases 
with | a; | from its zero value at x = 0. Also, the increase 
is more drastic with positive values of x. Indeed, only 
for x > Carnot efficiency is reached. For negative val- 
ues of x, rj max increases initially (up to x = —10 where 
?7max = 0.003) and decreases thereafter. 

Similar results were obtained by maximizing the effi- 
ciency at maximum power rj(uj max )- This is illustrated 



in Fig. [5] For x = 1, rj{u> max ) = rjc/2 and the Curzon- 
Ahlborn limit is recovered. From our results we find that 
at large asymmetry, f)(u) max ) decreases with the increase 
of |x|. It is interesting to remark that, in contrast to 
T)maxj the maximum value of T)((jj max ) is obtained when 
x 7^ 1. In particular, we find that very near to the sym- 
metric value x = 1, the Curzon-Ahlborn limit of effi- 
ciency at maximum power rj(u) max ) = rjc/2 can be ex- 
ceeded, as shown in the inset of Fig. [5] Even though 
there is only a small improvement with respect to the 
Curzon-Ahlborn limit, such result is interesting in that 
it provides the first evidence, in a concrete model, of the 
fact that the Curzon-Ahlborn limit, which is a univer- 
sal upper bound for time-reversal systems within linear 
response, can be exceeded when time reversibility is bro- 
ken. 

According to our numerical optimization, large asym- 
metry results in low thermoelectric efficiency. Also, we 
did not observe a significant improvement in efficiency by 
increasing the number of levels in the system up to six 
(for systems with a larger number of levels, it becomes 
difficult to obtain convergence). 



B. Gaussian Unitary Ensemble model 

To see whether thermoelectric efficiency can be im- 
proved in more general systems, we turn to Random 
Matrix Theory (RMT) models. Indeed, Hamiltonians 
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FIG. 3. Relation between maximum efficiency rj ma x and 
asymmetry parameter x [panel (a)] and between efficiency at 
maximum power rj{ui max ) and asymmetry parameter [panel 
(b)] for 20000 realizations of 3 x 3 Hamiltonian drawn from 
the Gaussian Unitary Ensemble. 



drawn from Gaussian Unitary Ensemble (GUE) give a 
good description of complex physical systems with bro- 
ken time reversal symmetry. To this end, in this section 
we study the efficiency (and asymmetry) distribution of 
GUE Hamiltonians. We consider a random Hamiltonian 
Hs drawn from a Gaussian Unitary Ensemble^. To com- 
pare with the results of the Aharonov-Bohm interferome- 
ter, we set the dimension of the matrix as n = 3. Results 
of one such calculation with 20000 different realizations 
are plotted in Fig. [3] Here, we take T = 0.1, /i = 0.1 
and 7L=7ii=7F = 0.1. Note that results are not sen- 
sitive to variations of the parameters T, /j,, jl, 7r, jp. 
Top panel represents the dependence of maximum effi- 
ciency rj max on the asymmetry parameter x and bottom 
panel the dependence of efficiency at maximum power 
f]( u max) on x. Results shown in the figure reveal that 
as the asymmetry in the thermopower is increased, the 
efficiency decreases. This implies that it is practically 
hard to achieve large efficiency with large asymmetry in 
thermopower and confirms our results obtained by opti- 
mizing the efficiency of Aharonov-Bohm interferometer 
in the previous subsection. 

Actually, we find that obtaining large values of effi- 
ciency is highly improbable. This is illustrated in Fig. [4] 
with the logarithmic distribution of i] max and rj(ui max ). 
Both r\ max and rj(uj max ) follows an exponentially decay- 
ing distribution indicating that large efficiency is expo- 
nentially improbable. On the other hand, there is no sig- 
nature of a cut-off or a sharp border, suggesting that in 
principle r) = f]c might be possible for some special choice 
of parameter values. Also, the distribution of asymme- 
try parameter x follows exponential distribution imply- 




1.50 



FIG. 4. Distribution of (a) maximum efficiency rjmax, (b) ef- 
ficiency at maximum power rj(u} max ) and (c) asymmetry pa- 
rameter x. Inset of (a) and (b) shows the distribution of 
rjmax and r)(cj max ), taking into account only data such that 
| a: — 1| > 0.5. Note that obtaining high values of both effi- 
ciency and asymmetry is highly improbable. 



ing that large asymmetry is also a rare situation in this 
model. Note that the maximum value of efficiency is ob- 
tained around the symmetric point x — 1. Hence, to 
better clarify the role of asymmetry we have plotted the 
logarithmic distribution of r] rnax and r)(u; max ), restricted 
to data for which \x — 1| > 0.5, in the inset of Fig. HJa) 
and HJb) respectively. From the figure, it is clear that 
distribution of r) max and r](Li max ) is exponential even in 
the absence of symmetry. 

In order to study the dependence of efficiency on the 
number of levels of the system, we studied GUE Hamil- 
tonians by increasing the dimension n. Our results with 
n = 4 and n = 5 (data not presented) show that there is 
little improvement in efficiency with the increase in the 
levels of the system. 

Thus, having studied both the Aharonov-Bohm inter- 
ferometer model and the broader class of GUE Hamil- 
tonians, we can conclude that for Hamiltonian models 
it is difficult to obtain large values of efficiency with 
large asymmetry. Hence, in the next subsection we turn 
to a more abstract transmission model and investigate 
whether the efficiency can be improved. 
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C. Transmission model 

Here, we assume that the information about the sys- 
tem is inaccessible and only the transmission probabil- 
ities Tij(E) at each energy E are known. Conservation 
of probability implies ry(I?) = 1, while the sum rule 
'^2 i Tij{E) — ^2jT~ij(E) ensures that the currents vanish 
at equilibrium. To simplify the calculations, the trans- 
mission probabilities are taken to be constant over a win- 
dow of energy. More precisely, each energy window is 
characterized by two parameters: Ek, the center of the 
window and A*,, the width of the window. Thus, in our 
model nj(E) = (cy) fc (0 < (cy) fc < 1) when energy 
E E [E k - A k /2,E k + A fe /2], and ry (E) = other- 
wise. Efficiency is optimized over the parameters tem- 
perature T, chemical potential fi, energies Ek, width A&, 
and transmission probabilities (cy) fc . We consider n non- 
overlapping transmission windows (k = 1, n). 

First, we consider a single energy window of transmis- 
sion (n = 1). Here, the integration domain R 3 in l|18p 
is restricted to the cube E x - Ai/2 < E { < E x + Ai/2 
(z = 1,2,3) and 



^12 — ^21 — A. 

x (E 3 
where the constant 



JJdEi 



df_ 

dEi 



E 1 ){E 2 -E 1 )(E 2 -E 3 ), (21) 



A = Sj {clp)i [(cpl)i + {cpr)i] {cpl)i • (22) 



The integral in (|2ip vanishes since it is an odd function 
of, for instance, E\ — E 3 and the integration domain is 
symmetric under exchange of E\ and E 3 . Thus, L' 12 = 
L' 21 and a symmetric thermopower is obtained. 

For n — 2 energy windows, we found numerically 
that the asymmetry parameter x is always limited to 
a finite interval. Infinitely large values of asymmetry 
are obtained for transmission with at least three en- 
ergy windows. For n = 3 windows, we need to opti- 
mize over 20 parameters: T, /j, Ek, Afc, and (cLR) k , 
(cLp)fc, (cpl),;, (cpfi)fc, (k = 1,2,3). (The other trans- 
missions (cij) k arc then determined from the conditions 

E i (cij)fc=E J -(c*j)jk = l;) _ 

Fig. shows the variation of the optimized maxi- 
mum efficiency r) max with the asymmetry parameter x. 
As in the case of the Aharonov-Bohm interferometer, 
Carnot efficiency r\c is obtained only in the symmetric 
case x = 1. Introducing asymmetry results in the reduc- 
tion of maximum efficiency from r\c- In particular, for 
positive values of asymmetry, r\ max increases quadrat- 
ically from zero value to rjc at x — 1 and decreases 
thereafter. When the asymmetry parameter is negative, 
Wmax increases with — x from x — to x — —50 (where 
Vmax — 0.24262) and then decreases for larger values of 
—x. 

Compared to the Aharonov-Bohm interferometer 
model, the rate of variation of maximum efficiency with 




X 



FIG. 5. Maximum efficiency r\ ma x versus asymmetry parame- 
ter x under three terminal transport with transmission model 
described in the text. Solid curve corresponds to the theoret- 
ically predicted upper bounds for rj ma x for a generic broken 
time reversal symmetry system. The optimized efficiency sat- 
urates the theoretical upper bound only near the symmetric 
value x — 1. 



asymmetry is much less for the transmission model. Also, 
the values of the optimized efficiency with large asymme- 
try are greater than in the Aharonov-Bohm interferom- 
eter model by more than one order of magnitude. For 
instance, when the asymmetry parameter x = 10, the 
optimized value of rj max = 0.004577c for the Aharonov- 
Bohm interferometer model and 0.2877(7 f° r the trans- 
mission model. To understand the large difference of 
efficiency obtained between the two models, we take a 
closer look on the optimized values of the parameters. 
Our analysis shows that optimal values of r] max are ob- 
tained when the transmission probabilities (cy), , which 
in principle could take any value between and 1, are 
either or 1. In particular, we find that the efficiency is 
maximal for particular combinations of the transmission 
probabilities, as described below. Out of the three com- 
binations (one for each window) , two correspond to sym- 
metric transmission {{cij) k — {cji) k for i.j — L,R,P) 
and one to perfectly asymmetric transmission ((cy). = 1 
and (cji) k — for some values of i and j). Also, for the 
symmetric cases only two reservoirs are involved in the 
transmission. For instance, the optimized values of pa- 
rameters for x = 1.4 are plotted in Fig. [6] Panel (a) 
corresponds to the derivative of Fermi distribution func- 
tion d *gjp with the optimized temperature T = 1 and 
chemical potential fi = 1.333. Three energy windows of 
transmission are 1) [-3.861, -3.713], 2) [-1.326,-1.325], 
and 3) [0.4931, 0.5201]. For the first case, t lr = t rl = 0, 
tlp = tpl = 0, and Tpp = Tp R = 1. Here, the trans- 
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FIG. 6. Optimized transmission functions for x = 1.4 using 
the transmission model. Panel (a) shows — /', with /' deriva- 
tive of the Fermi distribution function (f'(E) = ^* ), as a 
function of energy E with optimal temperature T = 1 and 
chemical potential /j, — 1.3333. The transmission probabili- 
ties tlr, trl, tlp, tpl, trp and tpr at each energy E are 
plotted in panels (b), (c), (d), (e), (f) and (g) respectively. 
Note that optimal transmission probabilities are either or 
1. 



mission is perfectly symmetric and is only between the 
right reservoir and the probe (i.e., R ^ P). Transmis- 
sion is perfectly asymmetric in the second energy window. 
This is clear from the figure as tlr = 0, trl = 1, t lp — 
l,Tpjj — 0, and trp — 0, trr = 1. Transport of electrons 
between the reservoirs is from L — s- R —> P — >• L. For the 
third energy window, the transmission probabilities are 
tlr = T RL = 0, T L p = T PL = 1, and t rp = t pr = and 
are completely symmetric between the left reservoir and 




FIG. 7. Efficiency at maximum 
metry parameter x under three 
transmission model described in 
dieted upper bound of rj(Lj max 
Dotted dashed line corresponds 
i)(w mal ) = Vc/2- Note that the 
ceeded in the interval [1,21 with 



power ^(umax) versus asym- 
terminal transport using the 
the text. Theoretically pre- 
) is plotted as solid curve, 
to the Curzon-Ahlborn limit 
Curzon-Ahlborn limit is ex- 
the transmission model. 



the probe (i.e., L P). We also obtained similar results 
using models with three delta peaks for transmission in- 
stead of transmission windows. Note that for the delta 
peaks model the width of transmission windows tends to 
zero i.e., Afc — > 0. 

Curzon-Ahlborn limit is slightly exceeded for small 
values of asymmetry in the Aharonov-Bohm interfer- 
ometer model. Can the efficiency at maximum power 
go significantly beyond the Curzon-Ahlborn limit? To 
check this we optimize efficiency at maximum power 
vi^max)- Results of our optimization are shown in Fig. 
LU As expected, Curzon-Ahlborn limit r](oj max ) — i]c/2 
is exceeded. When x > 0, r/(ui max ) increases quadrat- 
ically with increase in the asymmetry parameter until 
x = 1.349. At x = 1.349, it takes the maximum value 
of 0.56977c- Beyond 1.349, T}{u) ma x) decreases. More- 
over, for the asymmetry parameter x in the interval [1,2], 
vi^max) > Similarly to maximum efficiency, at 

large values of |a;|, the efficiency at maximum power de- 
creases with |sc|. Also, for x < 0, rj(u> max ) increases 
initially before the decrease (with the maximum value 
Vi^max) — 0.24261 obtained at x = —50). Note that 
for large values of asymmetry both the maximum effi- 
ciency and the efficiency at the maximum power are much 
smaller than the upper bounds set by thermodynamics 
(as discussed in Ref4) and shown as solid curves in the 
Figs. [5] and [7] 

We have also increased the number of transmission 
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windows, up to n — 7 (data not shown), and found only 
a slight increase of the maximum value of r](ui max ) up to 
0.570977c at n = 7. 

Our numerical results suggest that a three terminal 
thermoelectric transport yields large asymmetry with 
very low efficiency. Does the thermodynamics impose 
bound on maximum value of efficiency in non-interacting 
systems with broken time reversal symmetry under three 
terminal transport? Is it possible to achieve Carnot effi- 
ciency with asymmetry in the system under three termi- 
nal transport? These questions are addressed in the next 
subsection by considering a general model with random 
values for the Onsager coefficients. 

D. Random Onsager matrix model 



Even with all these restrictions, we find that a ran- 
dom Onsager matrix model can saturate the upper bound 
from thermodynamics for efficiencies (solid curves in 
Figs. [5] and [7]). In particular, rj(u) max ) — > r\c when 
|x| — > oo. However, it is clear from the definition of 
Onsager coefficients in Eqs. (fTS")) and (fTl)]) that for a non- 
interacting system there are further correlations between 
the different Onsager coefficients. For instance, the dif- 
ferent integrands are weighed by (E — fi) n ,n — 0, 1,2. 
For small systems, in particular for a three level system, 
these correlations are significant and put bounds on max- 
imum achievable efficiency as implied by our optimiza- 
tion results. We cannot exclude that more complicated 
transmission functions for a larger number of levels or 
interacting systems might allow approaching the Carnot 
efficiency even in the asymmetric case i^l. 



Consider then a 4 x 4 random Onsager matrix L. The 
elements Ly are chosen from the uniform distribution 
[—1,1]. In general, the elements of the Onsager ma- 
trix are not independent and are bound to fulfill certain 
conditions. First of all, it follows from thermodynamics 
that the entropy production rate is to be positive i.e., 
S = JiX 1 + J 2 X 2 + J 3 X 3 + J 4 X 4 > 0. This requires that 
the Onsager matrix L is positive-definite. 

Now, we turn our attention to specific case of non- 
interacting three-terminal system. In this case, there are 
further restrictions on the elements of Onsager matrix L. 
For instance, it is clear from the Eqs. (|T5|) and (TT5|) that 
-^13,^24,^31,^42 < (and that, as expected in general, 
all diagonal elements La > 0, i = 1, ...,4). Also, the off 
diagonal elements of each block matrix are even functions 
of magnetic field and hence are equal, i.e., Li 2 = £21, 
ii4 = L23, £32 = £41, £34 = £43- Note that this symme- 
try of the off diagonal elements is broken for the reduced 
Onsager matrix L', where we can have L' 12 ^ L' 21 . Posi- 
tivity of the transmission probabilities, Tji(E) > and of 
—df/dE, and the sum rule J2i T ij (E) = J2j T ij C®) imply 
£11 > I-Zasl) |£3i|; L 22 > \L 2i \,\L4 2 \; L 33 > |Li 3 |,|L 3 i|; 
and L 44 > |L 2 4|, 1^42 1 • 

There are further constraints on the values of Onsager 
matrix elements from the structure of reservoirs, which 
we assume to be ideal Fermi gases at temperature T and 
chemical potential [i. Upper bound on the energy inte- 
grals appearing in the Onsager coefficients are then given 
by 



eT 
~~h 

T 
h 



dE 
df(E) 



J dE (E-fx) rij (E) 
J dE (E - n) 2 Tij(E) 



T 



dE 
df(E) 



< 



< 



T 2 In 2 

T 3 7T 2 



24 



(23) 



The first and the third bound are saturated by setting 
Tij(E) = 1 for all values of energy, the second one by 
setting Tij = for E < \x and T\j = 1 otherwise, or vice 
versa Tu — for E > \x and tu = otherwise. 



IV. CONCLUSIONS 

We have investigated the thermoelectric efficiency of 
broken time reversal symmetry systems under three ter- 
minal transport. In these systems, the efficiency is de- 
termined by two parameters, the asymmetry parameter x 
and the "figure of merit" y. First, we have studied the op- 
timized efficiency of a realistic model of Aharonov-Bohm 
interferometer formed with three non-interacting quan- 
tum dots using simulated annealing. Our results show 
that Carnot efficiency r\c can be obtained when the ther- 
mopower is symmetric (x = 1). Introducing asymmetry 
in thermopower, the maximum efficiency r) max decreases 
from rjc- However, the efficiency at maximum power 
r)(cj max ) is maximal with asymmetry in thermopower 
{x 7^ 1). In particular, our studies illustrate that Curzon- 
Ahlborn limit can be exceeded in the linear response 
regime for a realistic model with broken time reversal 
symmetry. We have also studied the thermoelectric ef- 
ficiency of a generic model using random Hamiltonians 
drawn from GUE. Our analysis shows that it is highly 
improbable to obtain large values of both asymmetry in 
thermopower and efficiency. 

Furthermore, we have shown that the efficiency can be 
improved using an energy dependent transmission. In 
particular, optimizing the transmission matrix elements 
at three energy windows we have found more than one or- 
der of magnitude increase in the efficiency at large asym- 
metry over the three level Aharonov-Bohm interferom- 
eter model. The optimal values of transmission proba- 
bilities at each energy window are cither or 1. Also, 
two energy windows correspond to symmetric transmis- 
sion and one to perfectly asymmetric transmission. The 
Curzon-Ahlborn limit is exceeded in the interval [1,2] 
of the asymmetry parameter x and rj(uj max ) as large as 
0.57?7c is obtained. Similarly to the Hamiltonian models 
considered in this paper, the efficiency decreases at large 
values of asymmetry. 

Our extensive and accurate numerical results suggest 
that a three terminal thermoelectric transport is viable 
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only for large asymmetry with very low efficiency On the 
other hand, using a model with random values of Onsager 
coefficients one may obtain Carnot efficiency for maxi- 
mum efficiency T) max and efficiency at maximum power 
t)(w mM ) for arbitrary large values of asymmetry. How- 
ever, we argue that for non-interacting systems with a 
small number of levels there are correlations between the 
Onsager coefficients that bound the efficiency. In order 
to obtain large efficiency for large asymmetry, we need 



to turn our attention to interacting systems or to go be- 
yond the linear response. This remains to be analyzed in 
future. 

Note. After completion of our work, we became aware 
of a related work-22,, showing the existence of upper 
bounds on thermodynamic efficiencies for three-terminal 
transport as a consequence of the unitarity of the scat- 
tering matrix. The results from our transmission model, 
depicted in Figs. [5] and [JJ in practice saturate these new 
bounds. 
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